Saddles in the energy landscape: extensivity and thermodynamic formalism 



o 
o 

(N 



M. Scott Shell0 Pablo G. Debenedetti0 and Athanassios Z. Panagiotopoulotd 

Department of Chemical Engineering 
Princeton University, Princeton, NJ 08544 
(Dated: May 28, 2003) 

We formally extend the energy landscape approach for the thermodynamics of liquids to account 
for saddle points. By considering the extensive nature of macroscopic potential energies, we derive 
the scaling behavior of saddles with system size, as well as several approximations for the properties 
of low-order saddles (i.e., those with only a few unstable directions). We then cast the canonical 
partition function in a saddle-explicit form and develop, for the first time, a rigorous energy land- 
scape approach capable of reproducing trends observed in simulations, in particular the temperature 
dependence of the energy and fractional order of sampled saddles. 
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In recent years significant effort has been devoted to 
the study of supercooled liquids and their glasses Q] . An 
important aspect of these technologically significant 
systems is the interplay between dynamic and thermody- 
namic processes, which is thought to play a key role in 
their kinetic slowdown and eventual falling out of equi- 
librium at the glass transition 3J. The energy landscape 
formalism of Stillinger and Weber has been a useful tool 
in the theory of supercooled liquids Q . In this descrip- 
tion, a system's configuration space is partitioned into 
basins surrounding local energy minima. Termed "in- 
herent structures," these minima correspond to mechan- 
ically stable particle packings and are described statisti- 
cally by their depth: exp [Na(<p)]d<p gives the scaling of 
the number of distinct minima with per-particle poten- 
tial energy (or depth) <fc±d(/)/2, where N is the number of 
particles and a is called the basin enumeration function. 
Here, "distinct" refers to minima differing by more than 
mere particle permutation. This formalism permits a rig- 
orous transformation of the canonical partition function: 



D Nla(4,)-0<p-f3a vib (P,<j>)], 



(1) 



where [3 = 1/ksT and a V it>, the vibrational free energy, is 
the per-particle free energy when the system is confined 
to an average basin of depth (f>. For each temperature in 
the thermodynamic limit, the system samples basins of a 
well-defined energy </>*; deeper basins are accessed as the 
temperature decreases. One identifies a configurational 
entropy, NkBcr((f>*), which is that part of the entropy due 
to the multiplicity of amorphous configurations explored 
by the system. Good functionalities can be rationalized 
for a and a v ib (and measured in computer simulations) 
such that the partition function can be explicitly eval- 
uated this approach is also useful for characterizing 
kinetic processes as it has been observed that the config- 
urational entropy plays a key role in dynamics 0, ||| . 

Our work aims to extend the energy landscape formal- 
ism to include a description of higher order stationary 
points, i.e., saddles in the landscape Q. This approach 
provides a natural connection with dynamics 0, E| , but 



is more intricate than minima alone because, in addition 
to their energy, saddles are also classified by their or- 
der — the number of directions with negative curvature. 
In this Letter, we begin by deriving the extensivity prop- 
erties of saddles and propose their corresponding enumer- 
ation function. We then derive a saddle "equipartition" 
theorem and show the expected scaling behavior for low- 
order saddles. Finally we give the appropriate form of the 
partition function in this formalism and demonstrate its 
utility for describing the behavior of supercooled liquids. 

Our consideration of saddles relies on an extensive 
macroscopic potential energy. That is, a single macro- 
scopic system of N particles can be effectively divided 
into M <C N equivalent smaller subsystems, between 
which boundary interactions are negligible compared to 
the total energy. The number of particles in each sub- 
system is macroscopic, iVg ^> 1, but the number of sub- 
systems is also large, N$ <C N. This condition is satis- 
fied by most common types of molecular interaction (no- 
table exceptions include molecules with long-range inter- 
actions or which are themselves macroscopic in size) . For 
a single-component system of structureless particles, the 
potential energy can then be written as 



U(r N )^U(r^) + U(r^) + ... + U(r N s 



(M)> 
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where U is the potential energy function, r 
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r N } gives the positions of the particles, and 



{ r i,r : 

r^ s are the corresponding positions of the N$ particles 
in subsystem i. Here, any stationary point in the over- 
all system can be viewed as a combination of stationary 
points in each subsystem. It is relatively straightforward 
to determine the scaling behavior of minima (lfj: if the 
number of distinct minima in each subsystem is go, the 
total number due to their possible combinations is , 
or, exp [N In (go) /Ns] = exp [iVcoo] 11}. Our notation 
for the density-dependent constant a^o indicates its cor- 
respondence with the total number of inherent structures 
(not of a particular energy) ; equivalently, is the max- 
imum value of the basin enumeration function. 

For saddles, one must consider their order n, defined as 
the number of negative eigenvalues in the Hessian matrix, 
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, = Qr Qr ■ Imposition of Eq. gives rise to a Hessian 
which is reducible in each of the subsystems; therefore, 
the total number of negative eigenvalues is the sum of 
that for each of the subsystems. In other words, the sad- 
dle order in the total system is the sum total of the orders 
of the subsystems. We consider a particular distribution 
of the total saddle order n among the subsystems, letting 
the values Mi denote the number of subsystem saddles of 
order i = 0, 1 . . . dN$ (d = dimensionality). In this nota- 
tion, the constraints are J^i Mi — M and Yli * ' Mi = n. 
For an overall saddle order n and a particular distribution 
{Mi}, the number of distinct saddles is 



tt n ({Mi}) = (Ml/ JjAfj!) JJsf^ 



(3) 



where gi is the number of distinct saddles of order i 
in a subsystem. The total number of saddles, however, 
must be the sum of f2 n ({Mj}) for all possible distribu- 
tions {Mi}. We will find that one particular distribution, 
{Mi} max , overwhelmingly dominates this sum. First we 
switch to an "intensive" notation by introducing the fol- 
lowing variables: mi = Mi/M, and x = n/dN is the 
overall fractional saddle order. Insertion into Eq. [3] and 
application of Stirling's approximation yields 



n x (M,{m t }) 
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Similarly, the constraint equations in intensive form be- 
come Y2i m i ~ 1 ano - Si(*/^-^is) ' m i = x - These con- 
straints and the terms inside the brackets in Eq.^Jare all 
independent of the number of subsystems M. As a re- 
sult, the distribution {mj} max which maximizes the term 
in brackets depends only on the overall fractional saddle 
order x. In the thermodynamic limit, M — N/N$ — * oo, 
this maximum term dominates the sum over all distri- 
butions {mi}. These considerations lead directly to the 
saddle scaling behavior: 
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exp[7V6» 00 (x)] 



N/N s 



(5) 



where Cl x gives the number of distinct saddles of frac- 
tional order x. Here we have introduced the generalized 
function 8^ (x) which characterizes saddle scaling behav- 
ior and has the property 9 oc (x — > 0) = a^. Notice that 
the relevant order parameter for saddles is their fractional 
order, such that exp [N8 oa (x))dx gives the number of sad- 
dles with fractional order x ± dx/2. 

Returning to the subsystem scenario, we now find the 
distribution {mi} max which gives the dominant saddles in 
the system. To do so, one maximizes the term in brackets 
in Eq. 0] or equivalently, its logarithm. Using Stirling's 



approximation, the distribution {rni} max must satisfy 

max ^""^ mj In gj — m, In m, . (6) 

i 

By accounting for the constraints with the usual 
Lagrange multipliers and using Eq. [5] for gi — 
exp[Ns9oo(i/dNs)], the following form arises when Eq. 
|HJ is evaluated: 



exp [NsOM - jN s z] 
Ei' exp [NsOnW ~lN s z>] 



(7) 



where z = i/dN$ is the fractional order of a subsystem. 
Here 7 is a Lagrange multiplier ensuring the constraint 
Y]j z ■ mi = x. Because the terms in the exponential 
grow with the size of the subsystems, which are them- 
selves macroscopic, m, is essentially zero for all i except 
one, j max . With this simplification, the constraint yields 
*max = n/M, or z max = x. In other words, the total sad- 
dle order is distributed across the subsystems such that 
their fractional saddle order is equivalent to each other 
and to the overall fractional order. This is in effect an 
"equipartition" of saddle order across the geometry of 
the system. 

Such equipartition has important consequences for 
low-order saddles; it implies that the majority of such 
saddles are built from a collection of localized first-order 
saddles. In this sense, each direction of negative curva- 
ture in the potential energy corresponds to an elementary 
saddle "defect" in an inherent structure. This observa- 
tion can be used to determine the approximate behavior 
of Q x for small values of x. We assume that £dN non- 
interacting, first-order defects are possible for any inher- 
ent structure. The constant £, for example, may attain 
a value close to 1/d if each molecule in the system can 
move independently of the others so as to form a saddle. 
This yields 
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expiNcroo} x (£dN)\/(£dN - xdN)\(xdN)\ (8) 



or, taking the logarithm and applying Stirling's approx- 
imation, 

Mz^lWoo-f Inf - (l-f)ln(l-f). (9) 

We can also state for a system containing non-interacting 
first-order saddles that the average saddle energy de- 
pends linearly on order: 



$„ = $0 + n ($1 - $ ) 



(10) 



where is the average energy of an nth order saddle. 
This trend has been discussed previously in theoretical 
work and has been found in simulation studies to be 
appropriate [lil llfij. One must bear in mind, however, 
that both Ea.llUland the equipartition of saddles apply to 
the entire ensemble of stationary points, of which only a 
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minute fraction are sampled by the system in equilibrium 
at low temperature. A system may require, for example, 
the cooperative movement of many molecules in order 
to reach nearby saddles, which may cause low energy 
stationary points to be of higher order than expected. 
The low-T persistence of a linear relationship between 
saddle order and energy observed in simulation there- 
fore suggests that the first-order defect scenario persists 
even for low-energy saddles; cooperativity then arises be- 
cause the direction of negative potential energy curvature 
about these points is a superposition of several molecules' 
atomic coordinates. 

The considerations so far have categorized saddles only 
by their fractional order. Following the approach used for 
inherent structures Q, one might extend this descrip- 
tion to potential energy as an additional order parame- 
ter. We therefore introduce the saddle enumeration func- 
tion, 6{(j>i x), for which the expression exp [N9(<fi, x))d(j>dx 
is proportional to the number of saddles with poten- 
tial energy per-particle <f> ± d(f>/2 and of fractional order 
x ± dx/2. (The basin enumeration function is retrieved 
from 9{<p 1 x = 0).) This extension allows a meaning- 
ful casting of the canonical partition function in which 
configuration space is divided into "basins" surrounding 
saddle points 



o Nie(M-04-^ ib {M^) d< f, dx , (ii) 



In this equation, et v jb is the vibrational free energy around 
a stationary point of energy ip and fractional order x; 
formally it is given by 



f[U(r N )-Nj,] dj ,N 



(12) 
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where A is the thermal deBroglie wavelength, the aver- 
age is restricted to saddle points of energy (j> and order 
x, and the integral for a particular saddle k is performed 
over its associated configuration space Tk [l2j. In the 
large system limit, the integral in Eq. ^2 wm be domi- 
nated by the maximum value in the exponential, and the 
conditions for equilibrium can be written as 
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(13) 



The simultaneous solution to these equations provides 
the average saddle energy <f>* and order x* sampled 
by the system at specified temperature. The total 
Helmholtz free energy is then A/N = a V ib(A 4>* •> x *) + 
<t>* — kBT9((j)* , x*). One can identify from this equation a 
per-particle saddle entropy, ks9, which converges on the 
conventional Stillinger- Weber configurational entropy at 
very low temperatures when the system spends most of 
its time near minima and x* asymptotes to zero. This 
point of view sets the stage for a more rigorous connection 



with dynamics, for which x* contains pertinent informa- 
tion. 

Knowledge of the functions 9 and a v ib provides all the 
thermodynamic details of the system. We now show that 
a number of reasonable assumptions about their func- 
tional form results in a physically realistic and insightful 
picture. Our analysis addresses results for several well- 
studied glas s-forming systems for which numerical data 
exist pi. Il4. HH Hil Il7j . First, we assume the vibra- 
tional free energy to be independent of saddle energy, 
Gvib(/3, 4>i x ) w a v ib(/3, x). This is rigorously true for min- 
ima at absolute zero, but remains a working simplifica- 
tion in our analysis. Furthermore, we model the vibra- 
tional free energy in the classical harmonic approxima- 
tion 0: 
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(14) 



(3a vih « din (T^~ x Tg/T) - xlneA^/ f3 au l 2 



d\n(T s /T)-x[C[3-\n(irC[3)/2} (15) 



where as , cejj = \ are half the average curvatures 

of stable and unstable modes, respectively; Tq and Tjj 
are the corresponding Einstein temperatures [18j; I is a 
length scale characteristic of the saddle's associated con- 
figuration space volume [l9j |; and erfi is the imaginary 
error function. In the last line, we assume as = O-u 
and use the asymptotic expansion of the error function, 
erfi(a;) — * exp [x 2 ]/xy/ir, for the low T limit. C — ajjl is 
a lumped constant. In general, the harmonic approxima- 
tion is increasingly valid for stable modes at low temper- 
atures, but we have made liberal use of its application to 
the unstable modes, especially in that the final a v ;b de- 
pends nontrivially on the length scale I characterizing the 
e?A-dimensional integral. Nonetheless, this approxima- 
tion provides a starting point for analysis, and we leave 
investigation of more accurate forms to future work. 

For the saddle enumeration function, we assume a 
Gaussian form in energy, consistent with previous sim- 
ulation studies H [H H H2 : 



9(cb,x)= 9 00 (x) \ - {<(> - <j>{x)) /A 5 



(16) 



where (f> is the average energy of saddles of fractional 
order x and A is their characteristic energy range. For 
the dependence of the parameters in this expression on 
fractional saddle order, we use Eq. |^]for 9^, implement 
the linear relationship in Eq. ^] for <f> such that 4>{x) = 
4> + S ■ x, and assume A to be roughly constant. 

The usefulness of this approach can be seen in the pre- 
dictions of the theory. Such an analysis is possible by 
solving Eq. 1131 with the simplified expressions for the vi- 
brational free energy and the saddle enumeration func- 
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FIG. 1: Plot of the equilibrium fractional saddle order (x), 
and average saddle (Caddie) an d inherent structure (<f>is ) 
energies as a function of temperature. The squares, trian- 
gles, and diamonds are the respective results for a modi- 
fied Lennard- Jones system [l4| . Tmct ~ 0.435 is the mode- 
coupling temperature (Tui- 



tion. Wc choose representative parameters based on pre- 
vious simulations |l3| . though two have no precedent in 
the literature, which we instead choose so that x*(T) 
matches the simulation result in Rcf. Our results are 
shown in Fig. ^ Overall, the theory and assumptions 
capture the essential behavior observed in simulations 
(the general shape and relationship of the curves in Fig. 
^) and are a promising starting point for further refine- 
ment of the approach. We note that our parameters are 
chosen among several similar model systems, and quanti- 
tative agreement in Fig.E]might be possible if a complete 
dataset for any one system were available. 

In summary, we have presented a thermodynamic for- 
malism which includes higher order stationary points 
in the energy landscape. Through a reformulation of 
the canonical partition function and by using several 
physically-motivated simplifications, we show that this 
formalism captures important trends in the behavior of 
low-T glass-forming materials. Future work will investi- 
gate the relationship suggested by this approach between 
liquid kinetics and thermodynamics. 
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